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Abstract 

The BRAF (V600E) mutation in colorectal cancers that are microsatellite stable (MSS) confers a poor patient prognosis, 
whereas BRAF mutant microsatellite-unstable (IVISI) colorectal cancers have an excellent prognosis. BRAF wild type cancers 
are typically MSS and display chromosomal instability (CIN). CIN has not been extensively studied on a genome-wide basis 
in relation to BRAF mutational status in colorectal cancer. BRAF mutant/MSS [BRAFmut/MSS] cancers (n = 33) and BRAF 
mutant/MSI (BRAFmut/MS\) cancers (n = 30) were compared for presence of copy number aberrations (CNAs) indicative of 
CIN, with BRAFwM type/IVlSS (e/?/\Fwt/l\/lSS) cancers (n = 18) using lllumina CytoSNP-12 arrays. SR/AFmut/MSS and BRAFwXl 
MSS cancers showed comparable numbers of CNAs/cancer at 32.8 and 29.8 respectively. However, there were differences in 
patterns of CNA length between MSS cohorts, with BR/AFmut/MSS cancers having significantly greater proportions of focal 
CNAs compared to 6/?/\Fwt/MSS cancers (p<0.0001); whereas whole chromosomal arm CNAs were more common in 
e/?/\Fwt/MSS cancers {p<0.0001). This related to a reduced average CNA length in B/?/\Fmut/MSS compared to BRAFmMSS 
cancers (20.7 Mb vs 33.4 Mb;p<0.0001); and a smaller average percent of CIN affected genomes in 6R/AFmut/MSS compared 
to B/?/\Fwt/MSS cancers (23.9% vs 34.9% respectively). S/?/AFmut/MSI cancers were confirmed to have low CNA rates (5.4/ 
cancer) and minimal CIN-affected genomes (average of 4.5%) compared to MSS cohorts (p<0.0001). 6/?/\Fmut/MSS cancers 
had more frequent deletion CNAs compared to 6/?/AFwt/MSS cancers on 6p and 17q at loci not typically correlated with 
colorectal cancer, and greater amplification CNAs on 8q and 18q compared to BRAFwt/MSS cancers. These results indicate 
that comparable rates of CIN occur between MSS subgroups, however significant differences in their patterns of instability 
exist, with 6/?/AFmut/MSS cancers showing a 'focal pattern' and BR;AFwt/MSS cancers having a 'whole arm pattern' of CIN. 
This and the genomic loci more frequently affected in Bff/lFmut/MSS cancers provides further evidence of the biological 
distinctions of this important cancer subgroup. 
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Introduction 

The BRAF V600E mutation is present in approximately 10- 
15% of sporadic colorectal cancer (CRC) [1] and is a hallmark of 
the serrated neoplastic pathway of CRC, where cancers develop 
from serrated precursor polyps [2,3]. The CpG Island Methylator 
Phenotype (CIMP) is strongly associated with presence of the 
BRAF mutation [2,4,5]. In approximately half of these BRAF 
mutant cancers, CIMP related methylation and silencing of the 
DNA mismatch repair gene, MLHl, results in widespread 
frameshift mutations known as microsatellite instability (MSI). 
BRAF mutant/MSI cancers have been well characterized and 
show typical molecular and clinical features including an excellent 
patient outcome [4,6,7,8,9]. The remaining BRAF mutant cancers 



do not methylate MLHl and are microsatellite stable (MSS). These 
BRAF mutant/MSS cancers have not been as well studied, but 
importandy confer a very poor patient prognosis [9,10,1 1]. 

The majority of sporadic CRC are BRAF wild type and arise 
from conventional adenomas that follow a well defined pathway of 
molecular events leading to cancer [12]. These BRAF w'Ad type 
cancers are typically MSS and frequently show chromosomal 
instability (CIN) [8] , the presence of which has been correlated 
with a poor prognosis in these cancers [13,14,15,16]. Interestingly, 
the presence of the _Sit4_FV600E mutation in MSS cancers confers 
an even worse prognosis [9,17], however CIN has not been 
extensively studied on a genome-wide basis in this cancer 
subgroup. 
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CIN refers to the rate of acquisition of copy number aberrations 
(CNAs) where sections of DNA are affected by either deletion or 
amplification events [18]. CIN can affect whole chromosomes 
largely through dysfunctional chromosome segregation during 
mitosis [18,19], and aneuploidy is the stable state of abnormal 
chromosome numlx-rs [18]. Alternatively CIN can refer to the 
presence of widespread structural sub-chromosomal rearrange- 
ments resulting from incorrect repair of DNA damage [20]. These 
structural rearrangements can arise though repetitive rounds of 
breakage and fusion repair cycles leading to complex deletions, 
amplifications and translocations [21]. 

Few studies have extensively investigated CIN in the context of 
BRAF mutational and MSI status. We have previously found 
comparably high frequencies of LOH events between BRAF 
mutant/MSS cancers and BRAF wild type cancers at several key 
genomic loci (18q, 17p, 5q and 8p), that are known to harbour 
important tumour suppressor genes [22]. 

Apphcation of genome-wide single nucleotide polymorphism 
(SNP) arrays to study the presence of CIN has allowed the 
identification of different types of CNAs including complex 
aberrations and copy neutral loss of heterozygosity (cnLOH) 
events. Several common regions targeted by CNAs in CRC 
including deletions on chromosomes 17p, 18q, .5q, 8p, 4q and Ip, 
and amphtications on chromosomes 13q, 20q, 7p, 7q and 8q, have 
been confirmed through SNP array studies [23,24]. 

MSI and CIN have previously been considered as two distinct 
pathways of genomic instability due to findings of MSI cancers 
being largely diploid [13,25,26]. However, several studies using 
cytogenetic analysis have found MSI cell lines and cancers to have 
a considerable presence of chromosomal aberrations, predomi- 
nandy cnLOH events [27,28,29,30]. Similarly, studies have 
reported the presence of CIN and CIMP to be inversely correlated 
[31,32], and the incidence of frequent methylation to be associated 
with reduced rates and lengths of CNAs [33,34,35]. However, the 
majority of these studies did not stratify for presence of a BRAF 
mutation [31,33,35]. 

We and others have highlighted the importance of the 
BRAFmut/MSS cancer t\pe with their correlations with poor 
patient outcomes and presence of distinct molecular and clinical 
features [9,10,17,22,36]. This study expands on the characteriza- 
tion of these cancers by investigating the extent of CIN on a 
genome-wide basis which may help to determine further molecular 
aberrations that could be contributing to the aggressiveness of this 
cancer type. 

Materials and Methods 

Cancer Samples 

An initial cohort of 1052 sporadic colorectal cancers and 
matched normals were obtained from patients following surgery at 
the Royal Brisbane and Women's Hospital, Queensland, Aus- 
tralia. Written, informed consent was collected from all patients, 
and the study was approved under the RBWH and Bancroft 
Human Research Ethics Committee. Clinical data including 
patient gender, age, stage at diagnosis (American Joint Committee 
on Cancer, AJCC), and anatomical site of cancer (with proximal 
location considered as being proximal to the splenic flexure) was 
collected where available. 

BRAF, p53 and KRAS Mutation, MSI and CIMP Investigations: All 
cancer samples had previously been investigated for MSI status 
using the 5 marker panel of the National Cancer Institute 
(mononucleotide: BAT25, BAT26; dinucleotide: D5S346, 
D2S123, D17S250) and classified MSI if at least two markers, 
including at least one mononucleotide marker, were positive. The 



presence of the BRAF V600E mutation, p53 mutation (across 
exons 4 to 8), KRAS mutMion (at codons 12 and 13), and the CpG 
Island Methylator Phenotype (using a 5 marker panel consisting of 
CACAAIG, IGF2, jNEUROGl, RUNXS, SOCSl [5]) had also been 
previously determined [22,36,37]. The cancers were subserjuentiy 
divided into three cohorts depending on their MSI and BRAF 
mutational status: as BRAF mutant/MSS (n = 60), BRAF mutant/ 
MSI (n = 68) or BRAF wild type/MSS (n = 924). 

Single Nucleotide Polymorphism Arrays 

From these cohorts, 33 mutant/MSS, 30 _SA4_F mutant/ 

MSI and 18 BRAF wAd type/MSS cancers and matched normal 
samples were chosen for quantification using Picogreen dye, and 
anah'sis for genome-wide copy number aberrations (CNAs) with 
HumanC:)toSNP-I2v2.1 Single Nucleotide Polymorphism (SNP) 
arrays (lUumina; San Diego, Ca.) according to the manufacturer's 
instructions. The beadchips were scanned using lUumina's iScan 
system and the image data was analysed with lUumina's 
GenomeStudio version 2011.1.0.24550. The cancer traces were 
referenced to their matched normal profiles and the boundaries of 
all somatic copy number aberrations were manually determined 
and based on human genome build NCBI36/hgl9. To account 
for stromal contamination commonly present in cancer samples, 
tiie Simulated DNA Copy Number (SiDCoN) [38] and an 
automated SiDCoN2 tool which are R script based applications, 
were used to assign each CNA a log R ratio and B allele frequency 
score to determine the genotype of each CNA. The SiDCoN 
application was able to determine the extent of cells that showed 
an aberrant copy number change for each CNA, including 
heterogenous genotypes. Any individual CNA that scored less than 
20% of aberrant cellular involvement was excluded from analysis 
in order to ensure reliable CNA data [38,39]. The CNAs were 
then converted from Excel to custom data tracks and visualized on 
the University of California, Santa Cruz's Genome Browser 
(http://genome.ucsc.edu/) [40]. 

Cytogenetics terminology was appUed with 'gains' and 'losses' 
referring to whole chromosome arm events where a large genomic 
region was affected and typically consisted of small copy numbers. 
'Amplifications' and 'deletions' referred to CNAs covering sub- 
chromosomal or focal regions and these potentially involved 
greater copy numbers [41]. Specific types of deletion and 
amplification CNAs were analysed with deletion events comprising 
of loss of heterozygosity (LOH), copy neutral loss of heterozygosity 
(cnLOH) and homozogous deletion (HD) events, whilst amplifi- 
cation CNAs included 3n and ^4n (complex) amplification events. 

To identify the extent of CNA coverage per chromosome, the 
length of each CNA was calculated as a fraction of its coverage 
over the full length of the specific chromosome arm in order to 
allow for comparisons of CNAs occurring on all chromosome arms 
with differing lengths [41]. Continuous CNAs covering >95"/() of a 
chromosomal arm were termed 'whole chromosome arm' CNAs; 
regional CNAs covered between 50-94% of a chromosome arm; 
and focal events were considered as <50% the length of a 
chromosome arm in keeping with previously published data 
[24,41,42,43]. In this study, whole chromosome CNAs (continu- 
ous aberrations extending over both chromosome arms) were 
included in the analysis of whole chromosome arm CNAs as in 
Beroukhim et at [41] and the Cancer Genome Adas Network's 
characterization of CRC [24]. Minimal common regions (MCRs) 
were also identified and referred to the smallest genomic loci that 
contained deletion or amplification copy number changes at the 
highest frequencies across cancers in each cohort. 

Cancer Cell Density in Samples: SiDCoN assisted in estimating the 
cancer cell density of each sample [39], and those samples that 
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contained >40% of tumour cells were automatically included for 
analysis (Figure SI in File SI). As described by Dulak et al [43], 
the remaining cancers (5A4i? mutant/MSS 11/33 = 33%; BRAF 
mutant/MSI 12/30 = 40%; BRAF wild type/MSS 2/18=11%) 
were analysed for presence of co-existing molecular changes 
relating to tumourigenesis in order to confirm there was a 
sufficient ratio of tumour cell compared to normal cell content to 
justify molecular analyses. This analysis included the presence of 
methylated markers, evidence of MSI and LOH, and mutations of 
cancer-related genes [22,36] (Table SI in File SI). Data and 
statistical differences with these cancer samples either excluded or 
not were compared to verify their inclusion in this study (Table S2 
in File SI). 

Statistical Analysis 

Significant differences between categorical data were analysed 
with Pearson's chi-squared test, or Fisher's exact test where 
appropriate. Proportions were tested using a proportion test, and 
where appropriate these p-values were corrected for multiple 
comparisons using the Benjamini-Hochberg method. For contin- 
uous variables, ANOVA was used to test for a significant 
difference between groups, and Post-Hoc analysis (using Tukey's 
HSD) was performed to explore differences further. For tests 
within cohorts, either a paired t-test or Wilcox's sign rank test and 
Friedman's test of related samples was performed. P values £0.05 
were considered significant. 

Results 

Clinical and Molecular Features of Study Cohorts 

33 5A4^' mutant/MSS (BRAFmut/MSS), 18 BRAFwM type/ 
MSS (BRAFwt/MSS), and 30 BRAF mutant/MSI {BRAFmut/ 
MSI) cancers were analysed. The majority of BRAF mutant 
cancers derived from the proximal colon, whereas most BRAFwt/ 
MSS cancers were found distaUy (p<0.0001) (Table 1). The 
BRAFmut/MSS cancers presented mosdy at advanced stages 
(AJCC III and IV), compared to BRAFwt/MSS and BRAFmut/ 
MSI cancers (p = 0.03) (Table 1). BRAFmut/MSI cancers 
conferred a later age of onset compared to MSS cancers 
(p = 0.01. Molecularly, the CpG Island Phenotype (GIMP) was 
predominant in the BRAF mutant cohorts, particularly the 
BRAFmut/MSI cancers; whereas no BRAFwt/MSS cancers were 
GIMP high (p<0.0001) (Table 1). A7L4^ mutations were present in 



28% BRAFwt/MSS cancers and the mutual exclusivity of KRAS 
with BRAF mutations was confirmed. 

Rates of Copy Number Aberrations in Molecular 
Subgroups 

Individual copy number aberrations (GNAs) that had &20% 
cellular involvement as determined by SiDGoN [38,39] were 
included in the following analysis. Gancers that had less than 40% 
tumour content as estimated by SiDGoN [39] had substantial 
evidence of cancer related molecular changes (Table SI in File SI) 
[43], suggesting sufficient tumour cellularity to also detect GNAs. 
Statistical differences in the rate and type of GNAs occurring 
between cohorts remained valid when analyses were performed 
with their exclusion (Table S2). Therefore the full cohorts were 
considered for this investigation. 

The MSS cohorts had comparable average rates of GNAs per 
cancer with a rate of 32.8 per BRAFmut/MSS cancer and 29.8 per 
5A4^wt/MSS cancer. The BRAFmut/MSI cohort had a signif- 
icantly lower rate of 5.4 GNAs per cancer (p<0.0001) (Table 2). 

The average length of a single GNA in the BRAFwt/MSS 
cohort (33.4 Mb) was significantly longer than the average GNA 
length in the BRAFmut/MSS cohort (20.7 Mb) (p<0.0001), and 
the BRAFmut/MSI cohort (23.6 Mb) (p<0.0001) (Table 2). The 
length of each GNA was considered as a fraction of the length of 
the specific chromosome arm [41]. This showed significant 
differences in the average and median chromosome fraction 
affected by GNAs occurring between cohorts, with the BRAFwt/ 
MSS having the highest fraction of chromosome arm involvement 
compared to the BRAFmutant cohorts (p<0.0001) (Table 2). This 
difference in average GNA length corresponded to a greater 
average percentage of genomes affected by GNAs in the BRAFwt/ 
MSS cohort (34.9%; range 0-80.5%), compared to the BRAFmut/ 
MSS cohort (23.9''/o; range 0-68.6%). In comparison to MSS 
cancers, BRAFmut/MSI cancers had a minimal proportion of 
genome involvement (4.5%; range 0-25.8%) (p<0.0001) (Table 2). 
Due to this small extent of GNAs affecting BRAFmut/MSI 
cancers, the following results wiU mainly compare the two MSS 
cohorts. 

Deletion and amplification GNAs were considered as a rate of the 
total number of GNAs occurring within that cohort, as well as the 
number of events occurring per cancer within each cohort. Across 
all cohorts, deletion GNAs were more common than amplification 
GNAs, with deletion events constituting approximately 73% of all 



Table 1. Clinical and molecular data of cohorts. 







Feature 


e/?/1/iTiut/MSS 


BRAfVitlMSS 


BRAFmutmSl 


P Value 


n 


33 


18 




30 




Average Age at Onset 


68.5 


69.1 




76.2 


0.01 


Female Gender 


21/33 (64%) 


8/18 


(44%) 


22/30 (73%) 


0.13 


Proximal Location 


21/31 (70%) 


4/18 


(22%) 


27/29 (93%) 


<0.0001 


AJCC Stage 1 


1/25 (4%) 


3/18 


(17%) 


7/27 (26%) 


0.03 


AJCC Stage II 


9/25 (36%) 


8/18 


(44%) 


16/27 (59%) 




AJCC Stage III 


1 0/25 (40%) 


5/18 


(28%) 


2/27 (7%) 




AJCC Stage IV 


5/25 (20%) 


2/18 


(11%) 


2/27 (7%) 




CIMP High 


1 7/30 (57%) 


0/18 




21/30 (70%) 


<0.0001 


p53 Mutant 


1 2/30 (40%) 


9/18 


(50%) 


7/30 (23%) 


0.17 


KRAS Mutant 


0 


5/18 


(28%) 


0 





doi:1 0.1 371 /journal.pone.0091 739.t001 
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CNAs per cohort (Table 2). The BRAFwt/MSS cohort had a 
significantly greater average percentage of the genome afTected by 
amplification events than BBAFmut/MSS cancers (12.2% vs 5.2% 
respectively; p = 0.01) (Table 2). 

The most frequent deletion events occurring in at least 50% of 
cancers in both MSS cohorts, involved chromosomes Ip, 4q, 5q, 
17p, 18q and 22q. The BBAFmut/MSS cohort had significandy 
more common deletion events than the _BA4i'wt/MSS cohort at 
chromosomes 6p (p = 0.02), 6q (p<0.05) and 17q (p = 0.02) 
(Figure 1). Significantly more frequent amplification events 
occurred in BRAFwt/MSiS compared to BBAFmut/MSS cancers 
at 13q (p = 0.0009) and 7q (p = 0.006). The BBAFmut/MSS 
cancers had significantly more frequent amplification events at 8q 
compared to i?A4^wt/MSS cancers (p = 0.02) (Figure 1). 

The average number of the specific type of either amplification 
or deletion CNAs per cancer demonstrated the MSS cohorts had 
similar rates of types of events (Table 2). However the BRAFwt/ 
MSS cancers had significantly longer lengths of aU types of 
deletion and amplification events, except cnLOH CNAs (Figure 2). 
BBAFmut/MSI had significantly lower rates and shorter lengths of 
all types of events compared to the MSS cohorts (Table 2, 
Figure 2). 

Frequency of Copy Number Aberrations According to 
Length 

All CNAs were assessed for the fraction of coverage according to 
the specific chromosomal arm. Analysis of the length of all CNAs 
showed the vast majority were either less than 50% or longer than 
95 % the length of a chromosome arm for each of the three cohorts 
(Figure S2 in File SI). Therefore, in order to further compare 
frequencies of CNAs between cohorts, CNAs were considered as 
either whole arm (>95% chromosomal arm length), or focal (as < 
50% chromosomal arm length [24,41,42,43]). The remaining 
CNAs (50-94% chromosome arm length) were considered as 
regional events. Varying the threshold of the focal length from < 
35% to <65% chromosome arm length, still resulted in the 
majority of CNAs being kept in either the focal or whole length 
subsets, and did not alter the statistical significance of important 
findings (Tables S3A and S3B in File SI). 

Whole Chromosome Arm Copy Number 
Aberrations. The 5A4^wt/MSS cohort had a significandy 
higher propensity for whole chromosome arm CNA events at 32% 
compared to the 5A4i?mut/MSS cohort at 17% (p<0.0001) 
(Table 2). This corresponded to a significantly higher average rate 
of whole chromosome arm CNAs per cancer in BBAFwt/MSS 
compared to i^iMfinut/MSS cancers (p = 0.04); and a greater 
average proportion of genome affected by whole arm events in 
BRAFwt/MSS compared to BRAFmut/MSS cancers (22% vs 
12%, p = 0.02) (Table 2). The BBAFmut/MSI cohort had the 
lowest whole arm CNA rate of just 1 .4 per cancer (p<0.000 1); and 
just 3% of their genome affected by whole arm CNAs (p<0.0001) 
(Table 2). Within both MSS cohorts, the average number of whole 
arm losses were significantly greater than the average number of 
whole arm gains (BBAFmut/MSS p<0.0001, BBAFwt/MSS 
p = 0.001). BBAFwt/MSS cancers had significantly more whole 
arm gain events than BRAFmut/MSS cancers (p = 0.01) (Table 2). 

Regional Copy Number Aberrations. Rates of regional 
CNAs were siinilar between cohorts and whilst they occurred at a 
lower rate compared to whole arm and focal events, their inclusion 
allowed for a comprehensive description of CIN across all three 
cohorts (Table 2). Both MSS cohorts had significantly more 
regional deletion than amplification events per sample (Table 2). 

Focal Copy Number Aberrations. The BBAFmut/MSS 
cohort had the highest proportion of focal CNAs at 70.7% of all 
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Figure 1 . Percentage of cancers per cohort that had an amplification or deletion copy number aberration at each chromosome arm. 

Asterisl<s indicate those chromosome arms where significant differences (p<0.05) in the rate of CNAs per cancer occurred between MSS cohorts (red 
for the 6ft/AFmut/MSS cohort and green for the BRAFvjt/MSS cohort to indicate which has a significantly greater rate of CNAs per cancer). 
doi:10.1371/journal.pone.0091739.g001 



CNAs, whereas the BRAFwt/MSS had significantly less at 54.9% 
(p<0.0001). This equated to a rate of 23.2 focal CNAs per 
£A-Jftnut/MSS cancer, and 16.3 per £A4^wt/MSS cancer; the 
BRAFmut/MSl cancers had substantially fewer focal CNAs at 3.4 
per cancer (p<0.0001) (Table 2). There were significantly differing 
average lengths of focal aberrations per cohort with 6.3 Mb for 
BRAFmut/MSS, 10.9 Mb for BRAFwt/M&S and 2.3 Mb for 
BBAFmui/MSI cancers (p<0.0001). 

In aU cohorts focal deletion CNAs, predominantly through 
LOH events, were significantly greater than focal amplification 
CNAs {BBAFmut/MSS p = 0.004, BBAFwt/MSS p = 0.03) 
(Table 2), {BRAFmut/MSl p = 0.0001). Compared to BRAFwl/ 
MSS cancers, 5A4ftnut/MSS cancers showed significantly more 
frequent focal deletions at 18q (11/33, 33% Vs 1/18, 6%; 
p = 0.04) predominantly encompassing 18q2 1 .2 which includes the 
SMAD2 gene locus. Focal amplifications in iJA4iTiiut/MSS 
cancers were also more common compared to i?iJAFwt/MSS 
cancers at 8q (11/33, 33% Vs 1/18, 6%; p = 0.04) predominandy 



at 8q24.21 covering the Myc locus, and 18q (7/33, 21% Vs 0/18; 
p = 0.04) affecting 18qll.2 (containing G47:46' and CTAGE). 

Minimal Common Regions (MCRs). Minimal common 
regions were considered to include all lengths of CNAs. Both MSS 
cohorts showed a high rate of cancers (&40%) with targeted 
deletion events at several loci previously associated with CRC, 
such as 18q21.1-18q21.2 (which includes SMAD2, SMAD4, DCQ 
and 17pl3.1 (p53) (Table S4A in File SI). Loci not as commonly 
associated with CRC were also found to be deleted in a similar 
proportion of MSS cancers, and included 22ql2.1, 22qll.l, 
22ql3.2, 17pl2, 17pll.2, each of which contain several cancer 
related genes (Table S4A in File SI). 

Analysis of MCRs affecting &20% of cancers in at least one of 
the two MSS cohorts revealed several loci where the rates of CNAs 
differed substantially between cohorts. Although after adjustment 
for multiple comparisons significance was no longer reached, the 
5iL4ftnut/MSS cancers had a high frequency of deletion CNAs 
compared to BRAFwt/MSt^ cancers at several loci on 1 7q and 6p 




Figure 2. The average length of specific types of deletion and amplification copy number aberrations per cancer. There were 
significantly longer lengths for all events (except cnLOH) in the 6fi/AFwt/MSS compared to the B/?/\Fmut/MSS cohort. 6/?/lFmut/MSI cancers had 
significantly shorter lengths for all types of events compared to MSS cancers (p<0.0001). 
doi:1 0.1 371/journal.pone.0091 739.g002 
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BRAFmut/MSS BRAFwt/MSS BRAFmut/MSI BRAFmut/MSS BRAFwt/MSS BRAFmut/MSI 



Figure 3. Average number of copy number aberrations (CNAs) and percentage of genome affected per MSS cohort. A) Average 
number of CNAs delineated by length per cancer in each IVISS cohort. MSS cohorts had a similar number of overall CNAs occurring per cancer, 
however the 6/?/lFmut/MSS cancers showed a greater number of focal CNAs, with the BRAFwt/MSS cancers having a greater number of whole arm 
events. 6/?/\Fmut/MSI cancers had considerably fewer CNAs of all types. B) Average percentage of genome affected by CNAs delineated by length in 
each MSS cohort. BRAFwt/MSS cancers had the greatest proportion of genome affected by CNA events, which was due to the higher number of 
whole arm events in this cohort. 6R/Afmut/MSS cancers showed a lower proportion of the genome affected by CNAs, which is reflective of the 
comparably lower rate of whole arm and higher rate of focal events that occurred compared to BRAFwt/MSS cancers. 
doi:10.1371/journal.pone.0091739.g003 



including 17q22 (that contains cancer related genes RNF43 and 
VEZFl), 17q24.3 {S0X9) and 6p25.1 {CDYL) (Table 3). Amplifi- 
cation MCRs were more common in _6iMFmut/MSS than 
BRAFwt/MSS cancers at 8q24.21 [Myc), and 18qll.2 {GATA6, 
CTAGE) (Table 3). 

The BRAFmut/MSI cancers had substantially fewer MCRs 
than the MSS cohorts, however they did have a comparatively 
high proportion of cancers (&20%) with focal deletions at 3pl4.2 
{FHIT), 16pl3.3 {RBFOXl) and 20pl2.1 [MACR0D2) (Table S4C 
in File 1). 



Different Patterns of CIN exist between the BRAFmuXl 
MSS and e/?/AFwt/MSS Cancers 

Although the MSS cohorts had similar average numbers of 
CNAs per cancer (Figure 3A, Table 2), the 5AlPwt/MSS cancers 
had the greatest proportion of genome affected by CNAs 
(Figwe 3B). CIN in a typical £A4i<wt/MSS cancer predominantly 
occurred via whole chromosome arm events, whereas CIN in 
_BA4_Fmut/MSS cancers largely correlated with frequent focal 
CNAs which resulted in a smaller proportion of genome affected 
(Figures 3A and 3B). 




Figure 4. Heat map sKiowing the distribution of whole chromosome arm and focal copy number aberrations across the cohorts. 

Sample heterogeneity occurred within cohorts however a focal pattern is evident in the 6R/Afmut/MSS and a whole arm pattern is present in the 

BRAFwt/MSS cohort. 

doi:1 0.1 371/journal.pone.0091 739.g004 
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Figure 4 shows the genome wide distribution of CNAs across 
the three cohorts according to the type and length of event that 
occurred at a particular chromosome arm. Although the MSS 
cohorts do show sample heterogeneity, the predominantly focal 
pattern of CIN is evident in _SA4-Fmut/MSS cancers, and this 
contrasts to the whole chromosome arm pattern seen in BRAFwt/ 
MSS cancers. 

Discussion 

This study has shown that _BA4Fmut/MSS colorectal cancers 
predominantly harbour focal or targeted CNAs, whereas the 
BRAFwt/M&S colorectal cancers have significandy more frequent 
whole chromosome arm CNAs. This results in a greater average 
percentage of genome affected by CIN in fiAiPwt/MSS 
compared to _SA4Fmut/MSS cancers. BRAFwt/MSS cancers 
show a similarly high percentage of the genome affected by whole 
arm CNAs as a previous report of a large series of different cancer 
types, including CRC [41]. Comparatively, the _SAlFmut/MSS 
cohort has a significandy smaller proportion of their CIN affected 
genomes covered by whole arm events. Overall these observations 
identify that _SA4_Fmut/MSS cancers represent a more 'focal 
pattern' of CIN, whereas BRAFwt/MSS cancers display a 'whole 
chromosome arm' pattern of CIN. 

Across all cohorts, the frequency of deletion CNAs exceeded 
amplification events for all types of CNAs. This difiFerence may 
reflect a greater selection for deletions which could be tumour 
promoting and involve more simple mechanisms of acquisition, 
whereas amplifications may require more complex interactions 
with homologous and non-homologous chromosomes [44]. 

Whole chromosome arm CNAs were significandy more 
common in BRAFwt/MSS than in BRAFmut/MSS cancers. 
Whole chromosome arm CNAs can promote tumourigenesis 
through the gain of oncogenes and loss of tumour suppressors on a 
large scale [19]. However, whole chromosome arm CIN is also 
linked with cancer repression where the reverse of cancer 
promoting effects occur, and there is an overabundant loss of 
oncogenic factors and gain of tumour suppressive effects [19]. 
Additionally, increased chromosome copy number can lead to 
excessive protein production which may place greater metabolic 
stress on the cancer cell and ultimately reduce their rate of cellular 
growth and proliferation [19,45,46,47]. Wh(;th(;r the propensity of 
whole arm CNAs may be contributing to the less adverse nature of 
BRAFwt/MSS cancers compared to BRAFmut/MSS cancers 
through mechanisms described above, may warrant further 
investigation. 

The aggressive BRAFmut/MSS cancers had a significantly 
higher rate of focal CNAs across their genomes. There have been 
previous reports of early compared to late stage cancers 
harbouring more whole arm compared to focal CNAs [48], where 
stage I breast cancers were found to have more frequent whole 
chromosome arm CNAs compared to stage II/III breast cancers 
which had smaller, more complex events [49]. Furthermore, a 
detrimental clinical outcome in melanoma has been associated 
with a greater frequency of focal CNAS compared to whole 
chromosome arm events [50]. Potentially these complex, sub- 
chromosomal events may be facilitating cancer progression by 
specifically targeting key drivers of tumourigenesis. 

Different mechanisms relating to the origin of either whole 
chromosome or focal CNAs exist. It has been commonly reported 
that CIN involving whole chromosomes is due to errors relating to 
chromosome segregation during mitosis [19,51]. These errors are 
more likely to be those involving dysregulation of kinetochore- 
microtubule attachments, termed merotely, where a chromosome 



attached to both spindle poles mis-segregates at anaphase and 
results in whole chromosome aneuploidy [20]. The 'focal' pattern 
of CIN we have identified in BRAFmut/MSS cancers may 
associate with 'structural' CIN which involves structural sub- 
chromosomal rearrangements including deletions, amplifications 
and translocations [20] . The causes of these types of structural 
aberrations may involve dysfunctional repair processes of double 
strand breaks by homologous recombination and the error prone 
non-homologous end joining [52]. Potentially many of the 
particularly complex patterns of structural aberrations may not 
be drivc'r mechanisms in tumourigenesis but instead could be 
consequences of these disrupted DNA damage and repair 
processes. Studies of further genetic abnormalities unique to these 
specific CRC subgroups that could predispose to their respective 
patterns of CIN may be warranted. 

We have previously found that the CpG Island Methylator 
Phenotype (CIMP) and CIN can co-exist in BRAFmut/MSS 
cancers [22] . Potentially, the degree of methylation present and 
the subsequent effects on the extent of chromatin compaction may 
relate to the different rates of focal and whole arm CNAs observed 
between BRAFmut/MSS and 5A4Fwt/MSS cancers. Regional 
hypermethylation as present in CIMP positive cancers associates 
with increased levels of condensed chromatin, whereas widespread 
hypomethylation is present in cancers with a more open chromatin 
conformation [53]. This and our previous studies have found a 
substantial rate of BRAFmut/MSS cancers to be CIMP high 
[22,36], which may confer a more closed chromatin structure in 
these cancers. Global hypomethylation is well documented in 
CRC where it can associate with CIN and affects predominantiy 
BRAF wild type cancers [54]. A study found in regions with 
predominantly open or relaxed chromatin, repair mechanisms 
following double strand breaks were quicker to act due to greater 
accessibility of repair enzymes to the damaged site and 
subsequentiy resulted in less chromosome fragmentation [55]. 
These findings could help to account for the reduced rate of focal 
CNAs found in the BRAFwt/MSS cohort. The majority of 
BRAFmut/MSS cancers were CIMP high and demonstrated a 
'focal pattern' of CIN, which may suggest that a condensed 
chromatin structure contributes to a propensity of focal CNAs. 

CIN was evident in the majority of BRAFmut/MSl cancers, but 
affected a much smaller proportion of the genome compared to 
MSS cancers. Several genomic regions containing fragile sites, 
such as die FHITgene locus at 3pl4.2, RBFOXl at 16pl3.3 and 
MACR0D2 at 20pl2.1 that were targeted for deletion in MSS 
cancers, were also relatively commonly deleted in BRAFmut/MS\ 
cancers [56,57]. The lower degree of CIN present in MSI cancers 
may relate to findings that the onset of MSI is an early event in the 
development of MSI/CIMP positive cancers [58], and as this type 
of genomic instability is already present, there may be redundancy 
for the development of further genomic instability through GIN. 

As well as distinct variations in the pattern of GIN displayed 
between the two MSS cohorts, analysis of the minimal common 
regions (MGRs) of CIN revealed differential rates of either deletion 
or amplification CNAs occurring at certain genomic loci between 
them. Many of those more frequent in BRAFmut/MSS cancers, 
for example deletions at 6p25.l-6p21.33 and at specific loci on 
1 7q where several Wnt regulatory genes reside {RAF43, AXLN'2 
and S0X9) [59,60,61], have not commonly been associated with 
CRC. Additionally, BRAFmut/MSS cancers had a higher 
frequency of targeted amplification of 8q24. 13 that contains the 
Wnt signalling effector, Myc. These Wnt pathway related genes 
that may be specifically targeted in _Si?i4Fmut/MSS cancers could 
be an alternative mechanism in promotion of the Wnt signal in this 
cancer subtype. Amplification CNAs occurred at 18qll.2 in 
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BRAFmut/MSS cancers, whereas non-specific whole arm deletion 
events affected this region in BRAFwt/MSS cancers. The 18ql 1.2 
locus harbours two genes, GATA6 and CTAGE that have 
previously been reported to be amplified and upregulated in 
gastrointestinal cancers [39,43] including metastatic CRC [62]. 
Extended studies of these loci where greater rates of MCRs occur 
in _SiMi%nut/MSS compared to BRAFwt/MSS cancers may be 
indicated to ascertain whether these are driver mechanisms that 
may uniquely promote tumourigenesis in the BRAFmut/MSS 
cohort. 

This study has determined that a substantial presence of 
genome-wide CIN exists in the aggressive _SA4_Fmut/MSS cancers 
of the serrated neoplastic pathway. Significandy different patterns 
of CIN were found between the two MSS cohorts. BRAFnmt/ 
MSS cancers were found to harbour frequent focal length CNAs 
and therefore display a 'focal pattern' of CIN suggestive of 
commonly occurring structural rearrangements. Alternatively, the 
greater presence of whole arm CNAs in the BRAFwt/MSS cancers 
indicate they have a 'whole chromosome arm pattern' of CIN that 
may be due to dysfunctional mitotic events. Overall these findings 
suggest that either presence or absence of the BRAF V600E 
mutation could potentially affect subsequent acquisition of 
genomic instability in these subgroups of CRC. Extended studies 
to ascertain the clinical impact of the different patterns of CIN 
identified in these cancer subgroups may be warranted. Addition- 
ally, specific loci not as commonly associated with CRC that were 
more frequently affected by CIN in 5A4Fmut/MSS cancers were 
found, and this could help in the identification of molecular events 
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